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LOW FREQUENCY RADIO-FIR CORRELATION IN NORMAL GALAXIES AT ~ 1 KPC SCALES 
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ABSTRACT 

We study the radio-FIR correlation between the nonthermal (synchrotron) radio continuum emission 
at A90 cm (333 MHz) and the far infrared emission due to cool 20 K) dust at A70 /im in spatially 
resolved normal galaxies at scales of ~1 kpc. The slope of the radio-FIR correlation significantly differs 
between the arm and interarm regions. However, this change is not evident at a lower wavelength 
of A20 cm (1.4 GHz). We find the slope of the correlation in the arm to be 0.8 ± 0.12 and we use 
this to determine the coupling between equipartition magnetic field (Scq) and gas density (pgas) as 
Beq oc Pgas^^°'^^- This is close to what is predicted by MHD simulations of turbulent ISM, provided 
the same region produces both the radio and far infrared emission. We argue that at 1 kpc scales this 
condition is satisfied for radio emission at 1.4 GHz and may not be satisfied at 333 MHz. Change of 
slope observed in the interarm region could be caused by propagation of low energy (~ 1.5 GeV) and 
long lived (~ 10^ yr) cosmic ray electrons at 333 MHz. 

Subject headings: techniques: image processing - cosmic rays - dust - galaxies : ISM - galaxies : 
spiral - infrared : galaxies - radio continuum : galaxies 



1. INTRODUCTION 

The radio-far infrared (F IR) correlation in norma. 1 
galaxies was first observed by Ivan der KruitldTOTll [19731) 
and later extended by the IRAS mission. Subsequently 
it was established that the correlation holds good 
(within a factor of 2) over fi ve orders of rn a gnitude in 
radio and FIR luminosity (jCondonI 119921 : lYun et aD 
l2001t ) for a wide morphologic al class of galaxies like , 
spirals, irregu l ars and dwarfs (IWunderlich et ahl 119871 : 
iDressell 119881: iPrice fc Duri3 |1992D on global scales. 
Based on spatially resolved studies of normal and 
irregular galaxies it is seen that the correlation holds 
eve n at scales of few tens to hundreds of parsecs (see 
e.g IBeck fc GoUal fl98§. IXu et al .l [l99 2l: iHoernes et al.l 
1998t iMurgia et all 120051: iTab ataba ci et all l2007al: 



200l IGroves et all [200l . 
Cho fc Vishniad (I2OOO0 re- 



20O) 



Hughes et aLll2006l: iMurphv et al.ll20064 iPaladino et al 
20061: iPaladino. Murgia fc Orrul l2009t iDumas et al 



The basic model that connec ts these two regimes o f 
emission is via star formation ([Harwit fc Pacinil I1975D . 
The radio continuum emission arises due to synchrotron 
emission (henceforth nonthermal emission) from rel- 
ativistic electrons, produced in supernova remnants. 
A good fraction of them originate from massive (> 
10 Mq), short lived (< lO'* yr) stars. The FIR emis- 
sion arises from re-radiation by dust heated due to ul- 
tra violet (UV) photons emitted by the above pop- 
ulation of stars. Though the cause of the correla- 
tion is well understood, the tightness over several or- 
ders of magnitude still remains puzzling. Many mod- 
els explaining the correlation require close coupling be- 
tween the magnetic field {B) and the gas density (pgas) 
of the form, B oc P gas (see e.g., [^ lou fc Bi cavl 119931 : 
INiklas fc BecS [l997t iThompson et~l. 2006 ). Such a 
coupling can be established by magnetohydrodynamic 
(MHD) turbulence of the inters t ellar medium (ISM) (see 
IChandrasekhar fc Fermil [19531 ICho fc Vishniad [2000l: 



ICho. Lazarian fc Vishniad 
Numerical simulations by 
vealed that k = 0.5 is a manifestation of the equipartition 
condition, i.e, in steady MHD turbulence the magnetic 
field energy density and the energy density of the gas 
are similar. Similar values of k have been found through 
observations of magnetic field b y Zeeman splitti ng obser- 
vations in molecular clouds bv ICrutched (fl999( ). also by 
using equipartition magnetic field and molecular g as ob - 
servations in external galaxies bv INiklas fc Bec^ (|1997[) 
and in Milky Way and M31 bv iBerkhuiisenI (|1997l) . Al- 
tcrnatively, the slope of the radio-FIR correlation has 
been used to find k , where k ~ 0.4-0 .6 ("Nik las fc Beckl 
[1997: Hoc rnes et al.l 1998: Dumas et al . 2011). 

So far, spatially resolved and global study of the cor- 
relation has been done primarily using radio emission at 
1.4 GHz and higher freq uencies. The only low frequency 
study done at 150 MHz (jCox et al.lll988[) . confirms that 
on global scales the radio-FIR correlation holds good and 
is similar to what is seen at 1.4 GHz. To our knowledge, 
no low frequency (< 1.4 GHz, such as 333 MHz) spatially 
resolved study of the radio-FIR correlation exists in the 
literature. The motivation to do such a study arises from 
the fact that at lower frequencies the emission is largely 
nonthermal, hence better exhibiting the relation between 
magnetic field and star formation. Secondly, since the 
cosmic ray electrons (CRe) propagate larger distances in 
the galaxies at lower frequencies, it is important to assess 
how that affects the form of the radio-FIR correlation. 

In this paper, we present spatially resolved study of 
the radio-FIR correlation for four normal galaxies, NGC 
4736, NGC 5055, NGC 5236 and NGC 6946 at spatial 
resolution of ^1-1.5 kpc with radio observations made 
at 333 MHz (A90 cm) and 1.4 GHz (A20 cm). We also 
estimate the value of k and verify the equipartition as- 
sumptions. In Section[2]we discuss the various sources of 
maps used in this work and also define the parameter 'q' 
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TABLE 1 
The sample galaxies. 



Name Morphological Angular i Distance FIR Radio 

type size (D25)(') {°) (Mpc) A70^m A90cm A20cm 

(1) (3 (3) (4) (5) (6) (7) {&) 

NGC 4736 SAab 11.2x9.1 41 4.66^ SINGS GMRT Westerbork" SINGS (1374.5 MHz)^ 

NGC 5055 SAbc 12.6x7.2 59 9.2+ SINGS GMRT Westerbork SINGS (1696 MHz)* 

NGC 5236 SABc 11.2x11 24 4.51^ SINGS GMRT VLA*" GD array (1452 MHz)^ 

NGG 6946 SABcd 11.5x9.8 33 6.8=^ SINGS GMRT VLA C+D array (1465 MHz)^ 

In column (3) D25 refers to the optical diameter measured at the 25 magnitude arcsec"^ contour from lde Vaucoul curs et al] 1)19911 ). Golumn 
(4) gives the inclination angle (i) defined such that 0° is face-on. Distances in column (5) are taken from: ^ rKaracli cntsev et al.l (|2C!03), 
^ IKarachentsev et al.l 1I2OOI '). ^ IKarachentsev. Sharina fc Huchtmeieri 1 I2OOOI) and the NED t. Column (6) and (7) are the sources of data 
for the FIR maps and 333 MHz (A90cm) maps respectively. Column 8 are the data available at a higher frequency near 1 GHz (A20cm): 
IBraun et al.l II2007I) . ® VLA archival da ta using the CD array configuration (project code : AS325), ^ VLA archival map by combining 
interferometric data from C and D array, IBecM II2007I ) . 

"The Westerbork Synthesis Radio Telescope (WSRT) is operated 
by the Netherlands Foundation for Research in Astronomy (NFRA) 
with financial support from the Netherlands Organization for sci- 
entific research (NWO). 

'The Very Large Array (VLA) is operated by the NRAO. The 
NRAO is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 



which is used to quantify the correlation. In Section[3]we 
present our resuhs on spatiaUy resolved radio-FIR cor- 
relation using far infrared emission at A70 /im and radio 
emission at A20 cm and A90 cm. We discuss our results 
in Section m 



2. DATA ANALYSIS 

The fo ur galaxies in the sa mple for this study was cho- 
sen from iBasu et all (|2012[ ) . The large angular size of 
the galaxies ensure enough independent regions to carry 
out spatially resolved study. Our sample comprises of 
the galaxies NGC 4736, NGC 5055, NGC 5236 and NGC 
6946. Table [1] summarizes the salient features of our 
sample and the various sources of obtaining the archival 
data. 

To study the radio-FIR correlation using nonthermal 
radio emission, a thorough separation of thermal radio 
emission is needed. We used nonthermal radio continuum 
maps at A90 cm and A20 cm after separating the thermal 
free-free component mainly originating from Hii regions 
in recent star formation sites. Details of ob s ervati on and 
data analysis are discussed in IBasu et all ()2012|) . The 
thermal e mission was es timated using the technique de- 
veloped by iTabatabaei e t al. (2007b) , wherein, the dust 
extinction-corrected Ha map is used as a template for 
the thermal free-free emission. This is then extrapolated 
to the desired radio frequency and subtracted from the 
total emission map. The A90 cm maps were obtained us- 
ing the Giant Meterwave Radio Telescope (GMRT) while 
the A20 cm maps were obtained using archival data from 
various assorted sources (see Table [1]). The nonthermal 
maps had 40 arcsec resolution with 9 arcsec pixel size. 
We scaled the flux of each galaxy to a common frequency 
of 1.4 GHz using the spectral index map obtained from 
the 333 MHz and near 1 GHz images for each pixel. 

The galaxies were observed in the far-infrared by the 
Spitzer at A70 /zm as a part o f the Spitzer Infrared 
Nearby Galaxy Survey (SINGS; iKennicutt et"all 120031) 
using t he Multiband Im aging Photometer for Spitzer 
(MIPS; IRieke et al.l [200I . The images were obtained 
from the publicly available database in SINGS Data Re- 



lease 5^. The A70 ^im images have a pixel size of 4.5 arc- 
sec and a point spread function (PSF) of about 16 arcsec. 
These were convolved to the resolution of nonthermal ra- 
dio maps (40 arcsec) and re-gridded to a common pixel 
size of 9 arcsec. All the maps were then aligned to the 
same coordinate system. 

For the present study, the flux density per beam for 
the radio and FIR maps were determined within an area 
of ~ 40 arcsec diameter, with the adjacent region be- 
ing about one beam away to ensure independence. Pix- 
els with brightness above 2a {a is the rms noise in the 
map) were considered for the analysis. We estimate the 
slo pe of the radio-FIR a nd the quantity 'g' introduced 
bv lHclou. Soifer fc Rowa n- Robinson (198^. 

The parameter q is used as a measure of the radio-FIR 
correlation, where its dispersion indicates the tightness 
of the correlation. Conventionally it is deflned as the 
logarithm of the ratio of total FIR flux between A40 fim 
and A120 fim and the radi o flux measured at 1.4 GHz. 
However, we define q as per lAppleton et ahl (|2004[ ) using 
FIR flux density at A70 ^m, such that, 

qx = logio (Sro^nJSx) 

where, A is the radio wavelength (here, A = 20 cm or 90 
cm) and S'TOpm and S\ are the flux densities of A70 fim 
and radio wavelength respectively. The FIR emission 
from galaxies between A40 fim and A 120 fim is dominated 
by the emission fro m cool dust with dust temperature , 
Tdust ~20K (sec e.g.'Xu et al."1992VHoe rnes et al.|[T99l 
ITabatabaei et al . 2007b; Ba su et al. . 20121 ). The peak of 
this emission occurs at about AlOO fj,m. Note that a 
black body at ~20 K peaks at about A145 /xm, but a 
grey body [X~^ B\{T), where /3 = 2 is the dust emissiv- 
ity index and B\{T) is the Planck function) has a peak 
at ^ AlOO fiin. The maps at A70 ^m, which are nearest 
available to AlOO ^m, essentially traces this component 
of the dust. Using monochromatic A70 /zm emission to 
study the radio-FIR correlation does not affect the con- 
clusions sig niflcantly, except for a slight increase in the 
dispersion (|Murphv et al.|[2006aD . 

^ http: / / data.spitzer.caltech.edu/popular / sings / 
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Fig. 1. — Overlay of ~40 arcsec beams (marked in circles) on the Ha images" smoothed to 40 arcsec. The top and lower panels shows 
the arm and the interarm regions respectively (see Section 3 for details). 

"The im ages were download ed from the NED for the galaxies NGC 4736 (1-m Jacobus Kapteyn Telescope (JKT) at La Palma with filter 
Ha6570: [Knapen et al.|[200l ). NGC 5055 (2.3-m telescope at KPNO, filter: 6580) and NGC 5236 (0.9-m telescope at CTIO, fiher: 6563). 
For NGC 6946 the image was downloaded from the SINGS website. 



TABLE2 

Integrated flux densities of the galaxies at A90 cm (333 MHz: [Basu et al1|2012I) . A20 cm (1400 MHz; spectral index scaled 

FROM data given IN COLUMN 8 OF TABLE [TJ AND A70 flM (TAKEN FROM THE NED). THE MAP NOISE ((t)fOR THE 40 ARCSEC RESOLUTION 

IMAGES ARE ALSO GIVEN. 



Name 




f 90cm 


•S'20cm 


O" 20cm 


•S'70/jm 


O"70/Jm 




Jy 


mjy bcani"^ 


Jy 


mJy bcam~^ 


Jy 


mJy beam~^ 


NGC 4736 


0.9±0.06 


2 


0.31± 0.03 


0.35 


93.93± 7.34 


15 


NGC 5055 


2.3±0.13 


3 


0.41± 0.05 


0.4 


72.57± 5.16 


15 


NGC 5236 


6.86±0.62 


2.5 


2.36± 0.18 


0.3 


312.0± 15.6 


30 


NGC 6946 


4.3±0.24 


1 


1.5± 0.1 


0.2 


207.2± 16.1 


25 



TABLE 3 

Summary of the values of ij;^ and the fitted parameters for each of the galaxies. Columns 3, 4 and 5, 6 are the mean 

VALUES OF gx and THEIR DISPERSION AT A20CM AND A90CM RESPECTIVELY AS SHOWN IN FIGURE [2] THE VALUE OF q\ WAS COMPUTED 
USING THE FLUX DENSITY WITHIN ONE BEAM OF FWHM ~ 40 ARCSEC. COLUMNS 7 AND 9 ARE THE FITTED VALUES FOR IJA, WHILE 

Columns 8 and 10 are the slopes of the radio-FIR correlation. Here, a20cm, 620cm and agocm, 690cm are the parameters a 

AND b IN the fitted EQUATION, Sradio = ^.^IB ' WAVELENGTHS 20 CM AND 90 CM RESPECTIVELY. 



Name 




(g20cm) 


'''<320cm 


(lj90cm) 
























- Iog(a20cm) 


620cm 


- Iog(a90cm) 


690cm 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


NGC 4736 


arm: 


2.45 


0.12 


2.01 


0.20 


2.28±0.02 


0.73±0.04 


1.71±0.03 


0.49±0.04 




interarm: 


2.46 


0.15 


1.46 


0.23 


2.41±0.04 


1.13±0.12 


1.71±0.03 


0.28±0.04 


NGC 5055 


arm: 


2.26 


0.14 


1.59 


0.20 


2.14±0.02 


0.64±0.03 


1.44±0.02 


0.44±0.03 




interarm: 


2.03 


0.16 


1.27 


0.24 


2.10±0.01 


0.53±0.03 


1.38±0.01 


0.30±0.04 


NGC 5236 


arm: 


2.32 


0.14 


1.92 


0.20 


2.27±0.04 


0.92±0.04 


1.71±0.03 


0.66±0.04 




interarm: 


2.13 


0.20 


1.60 


0.26 


2.08±0.03 


0.76±0.05 


1.50±0.02 


0.40±0.04 


NGC 6946 


arm: 


2.31 


0.10 


1.81 


0.14 


2.26±0.02 


0.87±0.04 


1.69±0.02 


0.71±0.04 




interarm: 


1.97 


0.19 


1.32 


0.32 


2.07±0.01 


0.64±0.03 


1.52±0.01 


0.34±0.02 
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Fig. 2. — The figure shows the radio intensity vs. A70 fim FIR intensity (in Jy beam~^). The triangles arc for A90 cm and squares 
are for A20 cm. The filled symbols are for arms and unfilled symbols arc for interarms. The histograms are the distribution of ggocm 
and 520cm, where arms are shown with filled grey and interarms with unfilled histograms. The lines are the fit to the data of the form 
Sradio = o, X Si^Q^^ (See Table 3). The solid and dashed lines are fit to the arm and interarm regions. Black lines are for A90 cm, while 
grey lines are for A20 cm. 
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Fig. 3. — The percentage ehange in the value of {qx) between 
arms and interarms. Here, A(g) = (9A)arm — (QA)interarm for each 
of the sample galaxies. The triangles are for A90cm and squares 
are for A20cm. The mean change in the value of g20cm between 
arms and interarms is ~ 10 percent, while ggocm changes by ~ 30 
percent. 



3. RESULTS 

The spatially resolved study of the radio-FIR correla- 
tion was done by broadly classifying the emission from 
arm (including the central region) and interarm regions 
of these galaxies. The arms were identified from the Ha 
images for each galaxy. For the ringed galaxy NGC 4736, 
which has no prominent arms, the star forming ring was 
taken as the arm. The arm and the interarm regions used 
in our analysis are plotted as circles which are overlayed 
on the 40 arcsec Ha images in Figure [TJ The quantity 
q\ was computed within each such region. Note that the 
calibration uncertaint y at A70 fim could be ^ 20 per cent 
(|Murphv et al.|[2"006aD . This would lead to a systematic 
error of about 10 per cent in the values of qx. 

Table 2 gives the total flux density (in Jy) and map 
rms noise of the 40 arcsec images (in mJy bcam~^) of 
the galaxies. The galaxy integrated mean values of (720cm 
are 2.48±0.1, 2.25±0.07, 2.12±0.06 and 2.14±0.07, for 
the galaxies NGC 4736, NGC 5055, NGC 5236 and NGC 
6946 respectively. The qgocm are 2.02 ± 0.07, 1.5 ± 0.1, 
1.66 ± 0.09 and 1.68 ± 0.08 respectively. However, the 
spatially resolved estimates of (720cm and (^gocm suggests 
that their values vary between arm and interarm regions. 
Figure [2] shows the brightness of the nonthermal radio 
emission with the far infrared emission at A70 /im, both 
in units of Jy beam~^, for all the four galaxies. The figure 
also shows the distribution of q\ for A20 cm and A90 cm. 
It was seen that the star forming, gas rich spiral arms of 
the galaxies showed higher values for qx when compared 
to the adjacent low star forming interarm regions. Table 
3 summarizes the mean value of the quantity q\ and its 
dispersion for arm and interarm regions. 

The mean value in the arms for all the galaxies were 
found to be, ((720cm}arm = 2.32 with a narrow dis- 
persion of CTq^Q^j^ _arm = 0.14, while for the interarms 

(920cm) interarm = 2.15 and CTijjocm , interarm = 0.3. At 

A90cm, we find the ((?90cm)arm = 1-85 with cr, 



0.22 and (qgocm) interarm = 1-43 and (Jq 



g90cm,arm 



0.3. 



Figure [3] shows the percentage change in the value of 
qx between arms and interarms for each of the galaxies, 
where Aq = (qA)arm - ((7a) interarm- The squares are for 



A20 cm and triangles are for A90 cm. The mean of the 
'720cm changes slightly, by about 9 percent, between arms 
and interarms, however, this change is significant with > 
99.9 percent confidence using Kolmogorov-Smirnov test. 
At A90 cm, the change in the value of {qx) between arm 
and interarm regions are much pronounced with Aq ~ 
30 percent. 

3.1. Fit to the radio and IR flux densities 

The data were fitted using the form »S*radio — (2 x »5*jj^, 
where, ^radio is the flux density of the radio emission 
at A20cm and A90cm, and 5ir, is the flux density of 
the A70/xm infrared emission. The fitting parameter 'a' 



is an estimate for qx, where qx 



logj^Q a, such that 



logj^o Sx = ~q\ + b X log;^o Sir- The slope of the radio- 
infrared correlation is given by the parameter Sepa- 
rate fits were done for arm and interar m regions using or- 
dinary least-square 'bisector method' (jlsobe et al.lll990( ) 
in the log-log plane. The parameters obtained are sum- 
marized in Table 2 in the last four columns for A20 cm 
and A90 cm. The values of (720cm and (790cm are in good 
agreement with the mean values obtained from the re- 
spective distribution. The fitted parameters are plotted 
in Figure [H The black lines are for fits at A90cm and the 
grey lines arc at A20cm. The solid lines arc fit to the arm 
regions only and the dashed lines are for the interarms. 
All the correlations are highly significant in our case with 
Pearson's correlation coefficient, r > 0.8 (and r > 0.9 in 
most of the cases), except for interarm regions of NGC 
4736 at A90 cm, where r = 0.68. 

The slope of the A20 cm and A70 /im flux density for 
the arm regions for all the galaxies lies between ^ 0.65 
- 0.9. However, for the interarm regions the slope is 
slightly shallower, lying in the range 0.55 - 1. The mean 
value of the parameters for arm and interarm after the 
flt can be written as, 



logio S: 



10 'J20cm 



(2.24 ± 0.05) + (0.8 ± 0.08) logig ^ro^^m 
for arm (1) 



logio 520cm = -(2.17 ± 0.05) + (0.76 ± 0.14) logj^ ^yo^^m 

for interarm (2) 

The fltted value of (720cm and the slope differs slightly 
from arms to interarms. 

At A90cm, we find that the slope lies in the range ^ 
0.45 - 0.7 for the arms, whereas in the interarm region 
the slope lies in the range ^ 0.3 - 0.4. The slopes are 
much flatter than at A20cm. The mean values of the 
fltted parameters are found to be, 

logio 590cm = -(1.64 ± 0.05) + (0.6 ± 0.1) logip ^yo^m 

for arm (3) 



logio 5 



10 "-"OOcm 



^(1.53 ± 0.04) + (0.33 ± 0.07) logio ^yo^^m 
for interarm (4) 

There is a significant change in the value of (790cm and 
the slope between the arm and the interarm regions. 

3.2. q vs. ant 
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Fig. 4. — The figure shows the distribution of ctnt with g20cm (left) and (jgocm (right). The top panel shows the distribution of Ont 
estimated from A90cm and A20cm nonthermal emission radio maps within an area of 40x40 arcsec^ and a step size of 0.05, while the right 
side horizontal panel shows the distribution of qx within the same area and bin size of 0.05. The filled circles and histograms are for the 
arms and unfilled circles and histograms are for interarm regions for all the 4 galaxies combined. 



The quantity can be expressed in terms of CRe es- 
cape timescale (jesc) and synchrotron timescale (rsvn) as 
q\ (X log(rsyn/resc) (|Murgia et al.l 120051 iPaladino et al.l 
|2006[) . CRe emitting for these timescales also determines 
the variations in nonthermal spectral index^ (ant) from 
sites of injection to regions of dominant energy loss. In 
Fig. m we also study the variation of g20cm (left) and (jgocm 
(right) with the a^t determ ined at scales of ^ 1 — 1.5 kpc 
for all these galaxies from iBasu et ahl ()2012|) . The dis- 
tribution of a,it wit h steps of 0.05, ar e shown in the top 
panels (Figure 10 in IBasu et al.l[20T2l ). The right panels 
shows the distribution of qx for the respective frequen- 
cies with steps of 0.05. The shaded histograms repre- 
sents the arms and the unfilled histograms arc for the 
interarms. At A20 cm there is no apparent variation in 
the value of g20cm with ant (Pearson's correlation coef- 
ficient, r = 0.27), while at A90 cm the ggocm decreases 
with steepening of the ant (r = 0.7). 

In the arm regions where ant is seen to have a narrow 
distribution with mean —0.8 and dispersion of 0.14, the 
qx values are higher, suggesting Tgyn ^ Tosc. However, in 
the interarm regions, ant and qx have a wide distribution 
with more than 50 percent dispersion. The values of ^qx 
also systematically decreases as one moves from arms to 
interarms, indicating Tesc ^ Tgyn and thus the CRe loose 
energy before escaping the disk giving rise to steeper ant . 
Similar results were found for IC 342 a nd NGC 5194 
(jMurgia et al.l [20051 : IPaladino et al.l[2006l ). 

4. DISCUSSION 

We have studied the radio-FIR correlation at ^^1 kpc 
scales for four normal galaxies using nonthermal radio 
maps at A90 cm and A20 cm and the far infrared maps 
at A70 /xm. From the basic synchrotron theory (e.g., 
iMoffetl |1975() and considering the radio emission from 

^ In the text we define the spectral index a as, Sv ac 



CRe emitting at critical frequencies, the energy of CRe 
at A90 cm is ^1.5 GeV and at A20 cm is ^^3 GeV when 
they are gyrating in a typical magnetic field of ~ 10 fiG. 
The far infrared emission at A70 /im originates from cool 
dust at ~20 K heated by the interste har radiation field 
(ISRF) due to ~ 5 - 20 Up, stars (iDevereux fc EalesI 
[19891: [X3 [19901 : IXu fc Helod Il996t IDumas et al.ll20lH) . 
We separately examine these correlations for the arm 
and the interarm regions, that is, regions of high and 
low thermal fractions respectively. The results of the 
various parameters as discussed in Section |3] are given in 
Table 3 for individual galaxies, and here we discuss the 
average properties. The dispersion on the parameter qx 
is a measure of the tightness of the radio-FIR correla- 
tion, which for the arm region is found to be less than 
10 percent around the mean qx for both A20 cm and A90 
cm. For the interarm region the dispersion is seen to 
increase to around 20 percent for both the frequencies. 
Further we find the slope of the radio-FIR correlation 
for the arm regions (also the high thermal fraction re- 
gions) remains similar at both the radio frequencies (see 
Table 3). It should be noted that a large number of 
global scale radio-FIR correlation studies exist, where 
the observed slope is steeper and clos er to unity (see e.g., 
iPrice fc Duri3 119921 : lYun et al.l[200ll and the references 
therein). However, the spatially resolved studies relating 
FIR cool dust emission to A20 cm radio e mission, yields 
a val ue of the slope ~ 0.6 — .9 for LMC (Hughe s et ahl 
[2006[) and 0.80±0.09 for M31 ()Hoernes et al.a998iV It is 
difficult to compare the slopes obtained in global studies 
with the spatially resolved case. The flux in global stud- 
ies are averaged over both arm and interarm regions and 
we are uncertain about the contribution from each com- 
ponent. Multifrequency spatially resolved studies can 
provide an understanding of the relation between global 
scale and spatially resolved studies. For the present case, 
in the interarm regions (regions of low thermal fraction) 
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for A20 cm the slope is slightly flatter as compared to the 
arms (see Eq. 1 and 2). However, at A90 cm, the slopes 
become distinctly flatter than the arm regions (see Fig. [2] 
and Eq. 3 and 4). 

Our results can be used to determine the coupling be- 
tween magnetic field (B) and the gas density (pgas) as 
discussed in the introduction and thereby validating the 
'equip a.rtition' assumption s in these galaxies at 1 kpc 
scales. iDumas et al.l (|2011[ ) showed that the slope of the 
radio-FIR correlation relates to k as. 



n b 



3 - ant 
(n + l) b 

3 - ant 



optically thick dust (5) 
optically thin dust (6) 



where, n = 1.4ifc0.15 is the Kennicutt-Schmidt law index 
(see e.g., lKennicuttlll998() . b is the slope of the radio- 
FIR correlation and ant is the nonthermal spectral in- 
dex. For these face-on galaxies we use the assumption 
of optically thin dust to UV photons to estimate k. Wc 
find that k = 0.51 ± 0.1 at A20 cm and k = 0.4 ± 0.1 
at A90 cm. Similarly, for interarm regions due to a large 
range of ant we find k in the range 0.41 - 0.5 at A20 
cm and between 0.18 - 0.22 at A90 cm. Our estimated 
values of k, using the correlation between A20 cm and 
A70 ^m, are consistent with the predictions of numer- 
ical MHD simulations of different ISM tubulcncc mod- 



els, where k ~ 0.4 — 0.6 (see e.g.. Fiedler fc Mouschoviai 
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1993'; 'Ki m. Balsara fc Mac Lc^boOll: iThompson et al.l 
2006; .Groves et al.ll2003[) " 

In the arm regions, the slope and thus k remains sim- 
ilar for both A20 cm and A90 cm. Note that the above 
prescription to determine k is valid provided the radio 
and the FIR emission arises from the same emitting vol- 
ume, with a diameter of about 1 kpc for most of the 
observations reported here. In the arm regions the UV 
photon has a mean free path of ~100 pc within which 
most of the FIR emission arises. On the other hand, the 
CRe which gives rise to the radio emission diffuse farther 
away to - 1 kpc at 1400 MHz and - 2 kpc at 333 MHz 
in a galactic magnetic field of 10 /iG. Hence in order to 
have a similar slope with frequency, the energy spectrum 
of the CRe giving rise to the radio emission should be in- 
dependent of the volume element. This can only happen 
if the timescale for CRe diffusion/propagation (rdiff) is 
significantly larger than their generation timescale (Tgcn)- 
It turns out that the Tditr is about 8 x 10*" yr at 333 MHz 
and 4 x 10^ years at 1400 MHz which is significantly 
larger than the Tgen as evident from the supernova rates, 
which is one every 10^ — 10^ yr kpc~^ in Milky Way. We 
assume the same rate for these galaxies. 

The slope of the radio-FIR correlation in the interarm 
(low thermal fraction) region is similar to that of the arm 
at A20 cm, however it becomes distinctly flatter at A90 
cm. The flattening primarily happens due to relative in- 
crease in radio flux at A90 cm as compared to A20 cm, 
which has the effect that ant gradually becomes steeper 
in the interarms. This relative increase in A90 cm flux can 
be explained by continuous generation of CRe in the arm, 
which subsequently propagates into the interarm (e.g., 
from A to B or from farther regions in arms like C to B 
in Fig. O. The propagation timescale for these CRe are 
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Fig. 5.— The Ho image of the galaxy NGC 6946 (KPNO 2- 
m telescope, filter: KP1563) obtained from the ancillary data at 
SINGS website. The solid circles represents the diffusion scale of 
~1 kpc for ~ 3 GeV CRc at A20 cm, while the dashed circles 
represents the diffusion scales of ~2 kpc for ~1.5 GeV CRe at A90 
cm. See text for the details. 



few times 10^ years assuming Alfvcn velocity of 100 km 
and typical arm to interarm dist ance of 1-2 kpc . In 
such a scenario, using Equation 6 of (jKardas cv||1962|). in 
a typical galactic magnetic field of ^ 10 ^G, there would 
be a break in the energy spectrum for electrons above ~2 
GeV. This break frequency lies below A90 cm or above 
333 MHz. Such breaks have been seen at ~ 900 MHz 
and ~1 GHz for similar normal galaxies, NGC 3627 an d 
NGC 7331 respectively (jPaladino. Murgia fc Orr'^[2009( ) . 
Thus the CRe emitting at A90 cm, which lie above the 
break, do not loose significant amount of energy as com- 
pared to their higher energy counterparts. Hence, this 
results in increasing the relative flux at A90 cm. 

For the slope to remain similar between arms and inter- 
arm regions at A20 cm (below the break), the ratio of the 
radio to FIR flux densities should remain similar. Ob- 
served radio flux between arm and interarm changes by 
a factor of ~2-2.5. Similar ratio of flux density between 
arm and interarm regions at A20 cm can be caused due 
to steeping of the spectral index to < —1.1 as compared 
^ —0.6 to — 0.8 in the arms. This implies the FIR flux 
should change by a factor of ~^2.5-3 between arm and in- 
terarm regions for radio-FIR slope of ^0.8. The FIR flux 
density {F\) depends on the dust temperature (Tdust) 
and its density (pdust) as. Fx oc /9dustQabs(a, X)Bx{Tdust), 
where (5abs(«, A) is the FIR wavelength (A) d ependent 
absorption coefficient f or gain radius, a (jPrain c fc Le^ 
ll984HAlton et al.l[20M ). The temperature do not change 
sign ificantly between arm and interarm for these galax- 
ies (jBasu et al.ll2012|) . For a constant gas-to-dust ratio, 
i-e, Pdust oc Pgas, a factor of 2-4 drop in average gas den- 
sity between arm and interarm regions (found us ing the 
C0j:2-i.i maps from Heracles; iLerov et al.ll2009( ) would 
therefore cause the factor of 2-3 drop in FIR emission. 

The slope of 0.8±0.1 of the radio-FIR correlation indi- 
cates that the energy equipartition assumption between 
cosmic ray particles and magnetic field may be valid in 
the gas rich arms of the galaxies at our spatial resolution 
of ~1 kpc. For the interarm regions at A20 cm the slope 
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is similar to what is seen in arms, and thereby satisfying 
the equipartition conditions. The flattening of the slope 
at A90 cm does not indicate any break down of equipar- 
tition condition, but results due to overlapping emissions 
from adjacent regions. 

We thank Adam Leroy for kindly providing us the 
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Yogcsh Wadadekar for useful comments. We also thank 
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